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Abstract 

The fracture behavior of fiber-ceramics like C/C-SiC strongly depends on the initial dam¬ 
age arising during the production process. We study the transverse cracking of the 90° 
ply in [0/90]S' cross-ply laminates due to the thermochemical degradation of the matrix 
material during the carbonization process by means of a discrete element method. The 
crack morphology strongly depends on the fiber-matrix interface properties, the transverse 
ply thickness as well as on the carbonization process itself. To model the 90° ply a two- 
dimensional triangular lattice of springs is constructed where nodes of the lattice represent 
fibers. Springs with random breaking thresholds model the disordered matrix material and 
interfaces. The spring-lattice is coupled by interface springs to two rigid bars which capture 
the two 0° plies or adjacent sublaminates in the model. Molecular dynamics simulation is 
used to follow the time evolution of the model system. It was found that under gradual heat¬ 
ing of the specimen, after some distributed cracking, segmentation cracks occur in the 90° 
ply which then develop into a saturated state where the ply cannot support additional load. 
The dependence of the micro-structure of damage on the ply thickness and on the disorder 
in spring properties is also studied. Crack density and porosity of the system are monitored 
as a function of the temperature and compared to an analytic approach and experiments. 
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1 Introduction 


In the recent years extensive investigations were focused on enhancing high-tem- 
perature performance and reliability for fiber reinforced ceramic matrix compos¬ 
ites (CMC), especially their strength, damage tolerance and reliability as structural 
components. It was found that production process parameters control the forma¬ 
tion of the material’s micro-structure, and therefore, predetermine the activation 
of the micro-mechanisms, such as fiber-matrix debonding, fiber bridging and fiber 
pull-out, leading to damage tolerant behavior 0 - Therefore, one goal of ongoing 
research is to optimize the fiber-matrix interface properties in a way, that the inter¬ 
face is strong enough to avoid excessive fiber pull-out, but weak enough to activate 
fiber-matrix debonding, in order to avoid a catastrophic failure of the composite. 

The present paper focuses on carbon/carbon-silicon carbide (C/C-SiC). This mate¬ 
rial shows a distinctive crack pattern as a result of the production process, which 
can be influenced by the laminate layup as well as by the interface properties and 
the parameters of the production process itself. C/C-SiC is manufactured by the 
Liquid Silicon Infiltration (LSI) process [|^|^, that consists of three major steps. 
At first a carbon fiber reinforced plastic (CFRP) green compact is produced. This 
CFRP is then pyrolysed at temperatures ranging from normally 1100 K to 2000 
K, converting the thermosetting resin matrix into a glassy carbon one. In this step 
excessive cracking is observed, driven by the heavy shrinkage of the matrix, which 
is macroscopically prevented by the carbon fibers Q. Multiple defects like micro¬ 
cracks and micro-delaminations appear. To achieve good mechanical properties in 
combination with improved abrasive and oxidation resistance of the final CMC, 
the porous carbon/carbon (C/C) preform is infiltrated with liquid silicon at tem¬ 
peratures higher than 1684 K (melting temperature of silicon). Silicon reacts with 
both fibers and matrix to silicon carbide resulting in a nearly dense CMC with 
good overall properties [|^[^. The accessibility of crack systems for liquid silicon 
during the LSI process predetermines the later constitutive behavior of the C/C- 
SiC material in a dramatic way. A thermal fiber pretreatment can be employed to 
influence the interface shear strength and consequently the amount of degraded 
fibers [|6]-[^. In its extreme a siliconization of all fibers leads to a totally brittle be¬ 
havior of the composite. Therefore, the manipulation of the fiber siliconization is 
crucial for optimizing C/C-SiC with major influence on all engineering properties 
as well as on the prediction of the macroscopic failure behavior. Such efforts call 
for a comprehensive understanding for the development of the crack pattern during 
the pyrolysis and of the effects of internal structure on the macroscopic response. 
This paper presents an approach to understand the crack development during the 
pyrolysis of crossply-CFRP that combines a phenomenological study with a sim- 
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ulation based on a discrete element method. The phenomenological investigations 
are mainly based on thermo-microscopic observations of the crack development. 
Additional information was obtained by means of thermo-gravimetry, dilatometry 
and optical microscopy. 


Crossply material is built in thickness direction of alternating longitudinal (0°) and 
transverse (90°) plies. The mean topic of this study is the simulation of the fracture 
process during carbonization of crossply material and the characterization of the 
micro-structure of the occurring damage for different system sizes and disorder. 
Analytical and numerical predictive methods for the failure of the transverse ply 
were performed in the past at different length scales ranging from the microscopic 
scale, modeling single fibers [|9 ,10| and models for small fiber clusters pT| - [T3| 
up to the meso scale studying homogenized plies with a variety of model 
approaches, but mainly under mechanical loading. Small clusters of fibers with 
micro-structural disorder, embedded in a matrix material were modeled directly, 
employing fictitious or discrete crack [12| models. Unfortunately, these methods 
are very time consuming and results have limited meaning for the simulation of 
multiple cracking in thicker transverse plies. In addition, they do not allow predic¬ 
tions about size scaling and the possibility to model the dynamic fracture process 
itself. 


In the presented work, the 90° ply is represented by a two-dimensional triangular 
lattice of springs. The nodes of the lattice model fibers, oriented perpendicular to 
the plane of the lattice, and springs represent the matrix material in between. The 
spring-lattice is coupled by interface springs to two rigid bars which capture the 
two 0° plies or adjacent sub-laminates in the model. Disorder is introduced by as¬ 
signing randomly distributed breaking thresholds to the springs, i.e. a spring breaks 
if the load on it exceeds its breaking threshold. A detailed description of the model 
can be found in p3| |. To simulate cracking under thermal loading, internal load is 
imposed on the triangular lattice by varying the initial length of the lattice springs. 
The time evolution of the system is followed by solving the equation of motion 
of the nodes (molecular dynamics). With the discrete element method used in this 
study, realistic system sizes are possible, enabling simulation of the fracture pro¬ 
cesses with multiple defect interactions of thousands of cracks. Multiple cracking 
has been studied before, using random spring networks with springs of random 
distribution of strain to failure [16 -[M|, but rarely in the field of fiber composites 
research. This approach is particularly suited for studies on dynamic instability in 
crack propagation, the collective behavior of many interacting cracks and size ef¬ 
fects of multiple transverse cracking. Therefore, the discrete element method is 
practical for studying general features of the statics and dynamics of fracturing, 
like the crack morphology, global fracture patterns due to the collective behavior 
of many interacting cracks as well as the dynamic instability in the propagation of 
cracks. 


Our investigation is focused on the process of damaging, micro-structure of dam- 
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age, furthermore, on the relative importanee of damage meehanisms like segmen¬ 
tation and delamination, driven by the resin thermal degradation during the manu- 
faeturing proeess of the C/C from the CFRP, exposed to a moderate heating envi¬ 
ronment. Even though our model is eapable to prediet many aspeets of the fraeture 
progress (see (H), we foeus on miero-struetural properties that ean be eompared 
to experimental observations, like the development of segmentation eraek density 
and porosity. Sinee eonsiderable progress was made in understanding ehemieal 
kineties [T^20| and reaetion meehanisms [^22| of the earbonization proeess, 
this study is foeused on the ehallenge of eombining eomplex bulk material behav¬ 
ior, eomposite fabrieation and resulting global material behavior. Results from our 
model ean be used for optimization of the produetion proeess of all fiber eomposites 
that are pyrolysed as well as for the design of funetionally graded fiber eomposites. 


It is important to note that the gradual eraeking of thin film eomposites and eoatings 
exhibits eertain analogies to the gradual degradation of matrix materials during the 
produetion proeess of eomposites ( p3]-p^). In these systems a thin film or a so- 


ealled eoating is attaehed adhesively to a substrate and the sequential eraeking of 
the overlayer is studied either under meehanieal tensile loading or under thermal 
proeessing. Espeeially, the fragmentation of eoatings results in segmentation eraek 
seenarios whieh resemble very mueh to our system indieating that the underlying 
mieroseopie dynamies of eraeking is the same in the two oases. 


The paper is organized as follows: Seotion II gives an overview of the earbonization 
proeess, desoribing in-situ observations and measurements. The model oonstruo- 
tion, tests of the model and the numerioal results on eraeking obtained by simula¬ 
tions are presented in Seotion III. An analytio approaoh worked out in the literature 
to the thermal transverse eraeking problem is briefly summarized in Seetion IV and 
oonfronted with our simulation and experimental results. 


2 Experimental observations during pyrolysis 


The goal of the phenomenologioal researoh is to establish a hypothesis on the eraek 
development during the pyrolysis. This inoludes three major questions: (i) what 
types of oraoks do ooour? (ii) At what stage during pyrolysis do the different types 
of oraeks evolve? (in) How do the different types of eraeks interaet? Researeh was 
mainly done on laminates of fabrie, but for the study of speeifie questions, ex¬ 
periments were earried out on pure resin, uni-direetional laminates and eross-ply 
laminates. Sinee the simulation presented here is based on erossply material, this 
seetion foeuses on this type of reinforeement. 
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2. 1 Experimental procedure 


Experiments were performed on CFRP speeimens made from Tenax HTA fibers 
and a highly aromatie, nitrogen eontaining thermosetting resin by resin transfer 
molding. Uni direetional speeimens were produeed by winding the rovings on a 
steel eore before proeessing. Crossply material was made by using woven laminates 
with parts of the weft fibers being removed, reaehing a fiber volume fraetion of 
approximately 53%. All materials were eured at 473 K and post-eured at 510 K 
before proeessing them into samples. 


Thermo-gravimetrie (TGA) and dilatometrie analyses were performed to identify 
the main stages of the pyrolysis as well as to obtain the volumetrie ehanges. The 
porosity was measured by the water infiltration method [27 j after eooling the spee- 
imen from the given temperature. The eraek pattern after pyrolysis was observed 
with an optieal mieroseope with fluoreseensie eontrast. Note that only eraeks eon- 
tributing to open porosity are highlighted by this proeedure. In situ observations 
were performed under a mieroseope, equipped with the heating stage (see Ref. 0). 
The heating rate in all experiments was 10 K/min. 


2.2 Chemistry of the pyrolysis and its influence on mass and volumetric changes 


Chemieal reaetions start when the post-euring temperature at 510 K is passed. It ean 
be assumed that the post-euring temperature represents the stress-free state. This is 
used as a eriterion for the starting point of the simulation. Fig.[^ gives an insight in 
the degradation proeess by monitoring the mass losses together with the shrinkage 
of a pure resin sample. In general the pyrolysis ean be devided into four stages 
[ (28| : (i) During the post-euring stage up to 670 K mainly polyaddition reaetions 
oeeur with low mass losses and an expansion mainly ruled by the eharaeteristie 
thermal expansion of the matrix, (ii) Between 670 K and 920 K the main pyrolysis 
takes plaee, eonverting the polymerie matrix into a earbon one with dangling bonds 
and a deereasing amount of residual hydrogen, resulting in high mass losses, high 
shrinkages and an extremely high defeet density. (Hi) Up to 1470 K the remaining 
hydrogen is gradually removed with the observation of only low loss of mass (see 
Fig.[^). (iv) Beyond 1470 K mainly annealing reaetions and reaetions in the fiber 
oeeur [ [29| |, leading to further shrinkages in the eomposite speeimens (see Fig.[^). 
It is important to note that the glassy earbon from the pyrolysis of thermosetting 
resins is isotropie [(28 1. 


Fig. presents the ehanges in length of unidireetional and erossply samples in the 
laminate plane (Sx) and perpendieular to the fiber axis (efl. It is important to keep 
in mind, that all dilatometrie data are a superposition of the thermal expansion of 
the eonstituents and the thermoehemieal shrinkages of the matrix (eomp. Fig0. It 
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Figure 1. a) Thermo-gravimetric analysis together with dilatometer measurements on 
resin samples (mean value of nine measurements, three for each specimen orientation). 
b) Dilatometer measurements on UD and UD cross-ply specimen compared to numerically 


calculated strain for fiber volume fraction Vf = 0.6 from section 3.1.1 


was found that the stiff fibers result in expansion in fiber direetion that is one order 
of magnitude smaller than the out of plane shrinkage. Consequently we negleet the 
shrinkage in fiber direetion during the simulation. Note that the UD speeimen ean 
shrink in two dimensions, while the erossply one only has one direetion, leading 
to a higher shrinkage £z0°/90°. Additionally, a lower fiber volume fraetion of the 
erossply material eaused an enlarged shrinkage of 0°/90° in eomparison to the 
unidireetional 0°. More details ean be found in Q. 


2.3 Microscopic in situ observations of the pyrolysis 


The typieal miero-strueture developing during the pyrolysis is presented in Fig. 
Three major eraek types, fiber-matrix debonding (fissures without any preferred ori¬ 
entation), transverse eraeks and miero-delaminations are observed p0| . The eraek 
development starts with the onset of pyrolysis at approx. 670 K: evolving gases 
are trapped in the still eompaet CFRP, so they diffuse into pores leading to a net¬ 
work of fine ehannels (on the sub-mieron seale) through whieh the gases ean es- 
eape [22 311. The main eraek pattern is formed between 750 K and 900 K. First 
fiber-matrix debondings are observed followed by the eoaleseenee of debondings 
and pores forming transverse eraeks. The development of a eraek pattern formed by 
this eraek type is observed within a small temperature range. This ean be seen from 
Fig. [^, where the transverse eraek density of two thermo-mieroseopie speeimens 
is shown as a funetion of the pyrolysis temperature. Miero delaminations evolve 
from 800 K onward by eraek defleetion of the segmentation eraeks and lead to a 
saturation of transversal eraeking. As the general rule of eraeking during the pyrol¬ 
ysis, a linking meehanism was found: Debondings evolve by linking smaller eraeks 
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Figure 2. Side view on C/C cross-ply laminate after pyrolysis at 1173 K 

and pores while transverse eraeks and miero delaminations form debonded areas. 
This meehanism is in good agreement with the results of Sjoergren et al. [32| from 
investigations on meehanieally loaded erossply glass fiber reinforeed plasties. 


The development of the open porosity (see Fig. |^) refleets the development of the 
eraek pattern. The main inerease takes plaee during the main pyrolysis in the same 
temperature range in whieh the mesoseopie eraeks evolve. The maximum of 18% 
is reaehed after a pyrolysis up to 1170 K. Up to now, the meehanisms leading to 
the deerease of the porosity of 11% is not understood. 


3 The discrete element model 


For the study of the formation of defeets, a two dimensional triangular spring net¬ 
work model of eomposites eontaining rigid eells, is used, that is deseribed in detail 
in [ 151. Moleeular Dynamies (MD) simulation is used to follow the motion of eaeh 
eell by solving Newton’s equations of motion. In the present study we use a 4*^ or¬ 
der Gear Predietor Correetor seheme A general overview of MD simulations 
applied to eomposite materials ean be found in 


This method allows to simulate simultaneously a eonglomerate of eraeks within 
rather large lattiee sizes with even small eraeks being sharply defined. The funda¬ 
mental advantage of the lattiee model used in this investigation is its simplieity, 
giving direet aeeess and possible physieal interpretation to eaeh step of the algo¬ 
rithm. Consequently, one ean modify the rules of the model for features of interest, 
like the eharaeteristie properties of size, strength and of foree in a rather straight¬ 
forward and transparent way. 


The model is eomposed in three major steps, namely, (i) the implementation of the 
mieroseopie strueture of the solid, (ii) the determination of the eonstitutive behav¬ 
ior, and finally (Hi) the breaking of the solid. 


(i) Micro-structure: The model is eomposed of eireular eells of identieal radius 
Tf with their eenters loeated at the nodes of a regular triangular lattiee built of 
soseeles triangles of side length s, whieh is the eharaeteristie lattiee spaeing. The 
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Figure 3. a) Development of the crack density rho and b) the open porosity e', measured at 
room temperature. Since the in-situ observation of the development of p during the degra¬ 
dation process was limited due to low magnifications of the microscope, we scale these 
results to the saturation crack density evaluated for the whole specimen after pyrolysis. 

cells represent the eross seetions of fibers whieh are perpendieular to the model’s 
plane. The eell eenters are eonneeted by linear elastie springs, representing the 
material region with fiber, fiber matrix interfaee and embedding matrix material. In 
our simulation, the eells are the smallest partieles interaeting elastieally with eaeh 
other. We use a two-dimensional simulation with eells of unit thiekness, allowing 
only motion in the observation plane, with the two degrees of freedom being the two 
eoordinates of the position of the eell eenter of mass. The utilization of a regular 
lattiee is a elear negleet of the topologieal disorder, symptomatie for realistie fiber 
eomposite systems. Sinee we have no partieular basis for a topologieal disorder, we 
make a eomputational simplifieation by assuming the disorder in spring properties 
to aeeount for all relevant disorder present in the material. Disorder is given to the 
model by a Weibull distribution of breaking thresholds at the beginning of the 
simulation in the form 


P{Fd) = 1 - exp 



( 1 ) 


and is kept eonstant during the fraeture proeess (quenehed disorder). Fq represents 
the eharaeteristie strength of the material. The Weibull modulus m eontrols the 
degree of disorder in the distribution and it is usually ehosen in the range 2 < m < 
10, experimentally found to deseribe a variety of materials. The top and bottom row 
of eells, as well as the left and right rows are eonneeted to rigid beams with infinite 
stiffness via interfaee springs, representing the 0° — 90° interfaee of the eomposite. 

(ii) Constitutive behavior: The fibers are assumed to be rigid, undeformable eylin- 
drie bodies with the possibility to overlap when they have eontaet. This overlap 
represents to some extent the loeal deformation of the fibers (soft partiele approx- 
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Figure 4. Micro-structure of the model and cells with contact 

imation [351) and leads to a repulsive eontaet foree. The foree 'Pi on the eell i is 
the residual foree of all attraetive and repulsive forees on the eell ealeulated from 
their interaetion with neighboring eells. The eontaet law inside the system is only 
applied between two eircular partieles, after the breaking of the conneeting spring, 
thus preventing penetration of broken parts into eaeh other. It is a drawbaek of this 
method that the overlap force does not have a Hamiltonian formulation, which im¬ 
plies that no energy can be associated to the deformation modeled by the overlap. 
Stress is applied to the system by changing the initial length of the springs. The 
fiber elasticity modulus Eft, spring stiffness and breaking threshold is kept con¬ 
stant throughout the simulation. 


(Hi) Breaking of the solid: Stresses inside the system can be released either by the 
formation of surfaces (cracks) inside the transverse ply or by fracture of the inter¬ 
face elements on the boundary. In the framework of Discrete Element Methods, 
the complicated crack-crack interaction is naturally taken into account. If the total 
force on a spring F exceeds its breaking threshold (damage threshold), its stiff¬ 
ness is abruptly reduced to zero, resulting in a load redistribution to neighboring 
springs in the next iteration steps. If after some iterations the neighboring springs 
exceed their threshold value due to the additional load, they fail too, giving rise to 
crack growth. 


The time evolution of the particle assembly is followed by solving the equation of 
motion of the nodes (i.e. transverse fibers) 


mfi = Fi, i = (2) 

where N denotes the number of fibers of the transverse ply, and Ft is the total 
force acting on node T A 4*^ order Predictor-Corrector algorithm is used in the 
simulations to solve numerically the second order differential equation system Eq. 
Q. After each integration step the breaking condition is evaluated for all the intact 
springs. Those springs which fulfill this condition are removed from the further 
calculations (spring breakage). 
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3.1 Computer simulations 


Before applying the model to study transverse eraeking in fiber eomposites a variety 
of simulations have been performed in order to test the behavior of the model and 
to make parameter identifieation with respeet to experiments. Tests for the Young’s 
modulus, stress distribution between eraeks and Poisson’s ratio of the model system 
are deseribed in detail in [[T5|. 


3.1.1 Crack formation under gradual heating 


The effeet of gradual heating with the resulting shrinkage of the matrix from the 
degradation proeess is eaptured in the model by slowly redueing the initial length 
of the diserete spring elements, introdueing internal stress. The initial length of 
the spring elements is in the beginning the eharaeteristie lattiee spaeing s whieh is 
determined by the fiber radius and embedding matrix aeeording to Fig. as a 
funetion of the fiber volume fraetion vj 


j.f — 



Q<Vf< = 0.906. 


( 3 ) 


The initial length /°(T) is then ealeulated with the stress free thermal matrix strain 

£™(T) 

f{T)=2rf + {s-2rf){l-e^{T)). (4) 


During the simulation we linearly deerease the strain e^. Sinee the thermal strain¬ 
ing of the uni-direetional ply in fiber direetion is more than one order of magnitude 
smaller than in transverse direetion, we keep the x positions of the interfaee ele¬ 
ments fixed. 

The results are presented as a funetion of the dimensionless strain for the driving 
matrix shrinkage, giving general validity to the results. Later on we ean map the 
results for eomparative reasons to a temperature range of interest, using dilatometer 
measurements for the pure resin speeimen from Fig.[^, starting at approximately 
773 K. At this temperature the dilatometrie speeimens reaeh the length they had 
at the tempering temperature. Therefore, it is assumed that the thermal mismateh 
between longitudinal and transverse ply ehanges sign. The development of from 
the simulation shows good agreement with measurements in Fig. [^. At 773 K a 
regular network of pore ehannels has already developed. As the evolving gases 
ean now eseape the material, it is assumed that they do not have any influenee on 
the further eraek development and are therefore negleeted. TEM investigations on 
CFRP and C/C-SiC revealed a pore and pore ehannel diameter of less than 1/xm, 
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which is much less than the characteristie lattiee spacing s. The effect of the pore 
ehannels and small eracks and debondings of this size is therefore not taken into 
eonsideration explieitly, but through the disorder of the system expressed by the 
Weibull-statisties in Eq[^ 


The diserete element method gives the possibility to monitor the mieroseopie dam¬ 
age development in the speeimen. In aeeordanee with the experiments simulations 
show, that the damage proeess of the 90° ply is eomposed of two distinet parts. At 
the beginning of loading the weakest springs break in an uneorrelated fashion when 
they get over-stressed, generating miero-eraeks in the speeimen. This primary un¬ 
eorrelated micro-craek nucleation in Fig.[^ is dominated by the disorder of the sys¬ 
tem introduced for the damage thresholds of the springs. The micro-eraeks substan¬ 
tially change the local stress distribution in the ply leading to stress concentrations 
around failed regions, whieh gives rise to correlated craek growths under further 
loading. Growing craeks or craeks formed by the eoaleseence of growing eracks 
can span the entire thickness of the 90° ply resulting in segmentation, i.e. break-up 
of the ply into pieees, as it is illustrated in Fig. |^. Reaching the ply interfaee the 
eraek stops without the possibility of penetration into the 0° ply in the model, as 
it was observed in the experiments 0, but with the possibility of defleetion at the 
ply interfaee. Further segmentation eraeks mainly form between existing craeks 
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Figure 5. Snapshots of the model system of the horizontal and vertical size nx = 800 
(cells per row) and ny = 30 (number of rows), (a) nucleated micro-cracks distributed over 
the whole test cell, and the formation of the first segmentation crack. (6) Quasi periodic 
segmentation pattern, (c) Segmentation cracks and delaminations close to the crack density 
saturation state. 

until the crack density is saturated due to oeeurring miero-delaminations. Miero- 
delaminations usually start to oeeur when the eraek density has already reaehed 
high values, depending on the thickness of the transverse ply. Fig. shows that, 
as a result of the formation of spanning eracks, the ply gets segmented into pieces. 
The delamination zones along the interfaee can also be seen in Fig.[^. 

In order to obtain a quantitative eharaeterization of the miero-structure of damage, 
we analyzed the spatial distribution of miero-eraeks, i.e. the distribution of broken 
springs. A eraek is identified as a conneeted set of broken springs in the triangular 
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lattice taking into account solely nearest neighbor connections. A crack is consid¬ 
ered to be a segmentation crack if it spans from one side of the 90° ply to the other. 
Since cracks forming along the interface of the plies complicate the identification 
of segmentation cracks in the framework of our algorithm, first micro-cracks along 
the interfaces are left out of the analysis. Identified clusters with this method are 
demonstrated in Fig. [^for different values of disorder. The width of the clusters 
increases with the disorder, but still a good identification is possible even for large 
disorders. The method is illustrated in detail in [ [T5| . It can be observed that the 
number of micro-cracks and also the number of segmentation cracks increases with 
the strain and the segmentation cracks have a more or less quasi-periodic spacing. 
The formation of the segmentation cracks is very fast and often with stable crack 
growth only during the crack initiation. This behavior is in good agreement with 
the experimental observations. 



Figure 6. Snapshots of model systems of the size nx = 800, ny = 20 and identified 
clusters with disorder parameter m. 


The fraction of micro-cracks Ncracks/Nbonds and the dimensionless segmentation 
crack density Ngeg ■ t/l are presented in Fig. |^as a function of the strain je^ for 
several different values of the Weibull-modulus m. Sc is the average critical strain 
for spring elements. Since t is the ply thickness and I the overall system length, 
Ngeg 't/l can be understood as the thickness to width aspect ratio of segments. 
Varying m can be interpreted as one major effect of the thermal fiber pretreatment: 
decreasing m corresponds to a reduction of functional groups on the fiber surface, 
resulting in a clustering of functional groups along the fibers [361 and consequently 
in larger disorder in the system. The second major effect is the reduction of the 
interface shear strength as the number of chemical bonds in the fiber-matrix inter¬ 
face is reduced. It can be observed that cracking initiates at a finite strain value Sin 
called damage initiation strain, and starting from this point Ncracks monotonically 
increases during the entire loading process. Segmentation first occurs at a strain 
larger than after the number of micro-cracks reached a certain value. The num¬ 
ber of segmentation cracks Ng^g - t/l also increases, however, its value approaches 
a saturation value for larger strains in accordance with experimental observations. 
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Figure 7. (a) The total number of micro-cracks Ncracks divided by the total number of 
bonds Nhonds, and (6) the dimensionless crack density Ngeg ■ t/l as a function of strain 
jSc- The left row shows results for constant ply thickness of 20 cells and varying disor¬ 
der, while in the right row disorder is fixed to m = 4 and ply thickness is varied. Smooth 
curves were obtained by averaging over nine samples. 


Fig. and [^provide an insight into the role of disorder in the damage process. In¬ 
creasing the value of m, i.e. making the Weibull distribution narrower, the damage 
initiation strain Sin and the corresponding strain value for the start of segmentation 
increase, furthermore, the saturation value of the segmentation cracks increases. 
The slight decrease of Nseg for wider breaking threshold distributions is an artifact 
of the cluster analysis. A clear identification of segments becomes more difficult as 
segmentation cracks, located close to each other, are connected via delaminations 
and cracks. For large disorder it is possible that they are identified as one cluster, 
and consequently as one segmentation. 


3.1.2 Micro-cracking in laminates with varying cross ply thickness 

In order to get an insight into the dependence of the micro-structure of damage on 
the model thickness, systems were simulated with thickness ranging from approxi¬ 
mately 0.1mm up to 0.5mm. The micro-structure of damage in Fig.j^shows a quasi- 
periodic spacing for all calculated thicknesses. For the thicker models branching of 
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Figure 8. Micro-structure of damage for varying model thickness. The number of rows of 
cells ny is given in the figure. 

the segmentation craeks is observed elose to the ply interface. This effect can com¬ 
plicate the cluster evaluation, since neighboring segmentation cracks can coalesce. 
Thinner systems reach their saturation crack densities at higher strains than the 
thick ones. These results are compared to an analytical approach and experiments 
in section m 


3.1.3 Development of the porosity 

The porosity of the pyrolysed material governs the constitutive behavior of the C/C- 
SiC material, since the infiltration with liquid silicon mainly depends on the acces¬ 
sibility of the crack systems from the outside. This part of the porosity is called 
open porosity, while the total porosity is defined as the total volume of cracks, 
including those located deep inside the material without any connection to the out¬ 
side. Therefore, all elementary cracks that are not connected to the outside have to 
be identified in a cluster analysis. These isolated cracks are used to calculate the 
area of the closed porosity by adding the crack area of all broken elements 
forming theses cracks. The crack area between element i and j is approximated 
with equilateral triangles, computing the difference between the actual area and the 
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uncracked, stress free area at this temperature, 




( 5 ) 


Since segments delaminate from the ply interface in the course of the simulation, 
the open porosity can be evaluated by using the difference of the areas 

Aoc{T) = A{T) - A\T) - A^^T), (6) 


with A^(T) being the reference value of the stress free system without any damage 
and consequently with a porosity of 0%. A^^ is the sum of the area of all broken 
elements ij of all identified isolated cracks. The total area of the model A{T) is cal¬ 
culated from the positions of all edge elements. It is obvious that the open porosity 
calculated before stress can be released by forming cracks is artificial, but as soon 
as crack patterns form, the assumptions for the calculations become realistic. Nev¬ 
ertheless, the open porosity from measurements and from optical area analysis (see 
Fig. 1^ is approximately 13 % and therefore larger than the calculated ones that 
are 6-7% for this case. Even though quantitative agreement is beyond the scope of 
the model, the effects of disorder and ply thickness on the porosity are obtained 
correctly. Note that the two-dimensional model, which represents a small zone in¬ 
side the material, is limited to its characteristic scale and can not take effects into 
account, that are caused by the third dimension. A pyrolysed resin sample as well 
as a fiber bundle, therefore, are calculated with the porosity equal to 0%. On the 
other hand, the real material is irregular with resin rich zones, small delamination 
zones and a large specimen surface. It is also possible, that water intrudes the small 
cracks and gas channels along the fibers below the model scale. The dependence 
of the total porosity on the disorder of the system in Fig. is very small, but for 
the system with low disorder, the closed porosity (Fig. |^) is smaller than for the 
highly disordered system. Even though thick systems have a much smaller number 
of segmentations cracks than thinner systems, the open porosity in Fig. ^ does 
not show a strong dependence on the system thickness. This is mainly due to in¬ 
creasing crack openings in thicker systems. The closed porosity shows a stronger 
dependence on the thickness, with increasing values for increasing system sizes. 
It is therefore responsible for the deviation of the total porosity in Fig. ^ from a 
master curve. 


4 Analytical model 


The micro-mechanics of damage (MMD) analysis, used here for comparative rea¬ 
sons, is described in detail in p7| and was employed in [15|. A MMD analysis 
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Figure 9. Development of the porosity for varying disorder and eonstant model thiekness 
ny = 20 (a — c) and varying model thiekness but eonstant disorder m = 4 (d— ) respee- 
tively. 

is based on the individual prediction of the initiation and propagation of the vari¬ 
ous types of damage, including possible interactions between all identified relevant 
damage mechanisms. The analysis consists of two steps, (i) a stress analysis, and 
(ii) a damage analysis for the expected damage mechanisms. The stress analysis 
has to be undertaken in the presence of observed damage and is completely inde¬ 
pendent from the postulate of a failure criterion for the initiation and evolution of 
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damage, but of course not from the damage itself. 


(i) Stress analysis: The stress field in the presence of micro-cracking and delami¬ 
nation damage is calculated for a unit cell bounded by cracks located at each end 
(Fig. [T0| ). The entire specimen is built by a sequential formation of unit cells. A 
two-dimensional analysis for the x — z—plane derived in Ref. p8| using a varia¬ 
tional mechanics analysis is employed. Following [381 we assume that the x-axis 
tensile stress in each ply group is independent of 2 ;, that cracks span the entire ply 
thickness and that delaminations are always symmetric to the crack plane [[T7||. The 



Figure 10. Side view on a unit cell representing a cross ply laminate with segmentation 
cracks (a=half crack spacing) and delaminations di , d 2 in the transverse ly with half thick¬ 
ness ti. 

x-axis tensile stress of a micro-cracked laminate can be derived with force balance 
- 'ipix) (7) 


with the undetermined function the applied global stress cxo and A = ^ 2 /^ 1 , 
expressing stress perturbations caused by the segmentation cracks. The integration 
of the stress equilibrium equations with the unit cell boundary conditions gives 
the shear and transverse stresses in the unit cell in terms of ^jJ{x). The function 
^p{x) that minimizes the complementary energy gives the best approximation to the 
micro-cracked cross-ply laminate stress state. The solution of the Euler equation 
for finding ^p{x) is dependent on the temperature and can be found in detail in [381. 


(ii) Damage analysis: The failure criterion used, is an energy criterion, allowing 
a formation of a new segmentation crack, whenever the strain energy release rate 
(ERR) Gm associated to the formation of a new micro-crack exceeds some critical 
value Gmc^ also called intra-laminar fracture toughness or micro-cracking fracture 
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toughness. The ERR is also ealeulated with the two-dimensional, variational me- 
ehanies analysis [371. Miero-eraek indueed delamination eompetes with the forma¬ 
tion of new miero-eraeks. Experimental and numerieal studies show, that miero- 
eraek indueed delamination does not start before a eritieal miero-eraek density is 
reaehed. 


The two-dimensional, variational meehanies analysis deseribed above ean be ex¬ 
tended to aeeount for delaminations emanating from miero-eraek tips. The assump¬ 
tion that is independent of 2 ;, indieates identieal delamination lengths emanat¬ 
ing from the top and bottom miero-eraek tips. Eollowing [37| we furthermore as¬ 
sume that Gmc and the eritieal ERR assoeiated with the growth of delamination 
Gdc are equal. This way, the dimensionless eraek saturation density in Eig.[TT]was 
ealeulated. In principal, it is possible to get better agreement between MMD anal¬ 
ysis and the experimental data by adapting Gmc and Gdc until a satisfying fit is 
achieved, since these values have an ample scope for physical interpretation. Nev¬ 
ertheless, we find qualitative agreement in the increase of f/a as a function to the 
90° ply thickness for numerical and analytical results (compare Eigs.[^and |^. 


Thermal ply degradation with MMD: During the carbonization process, the shrink¬ 
age is opposed to the thermal expansion of the material. Therefore, experimental 
data on the strain-temperature behavior of the material from Eig.l^is required. In- 
situ experiments on the development of material properties like Young’s modulus 
Ex,Ez or material strength during the carbonization process without influencing 
the process of crack formation is believed to be beyond todays possibilities. We use 
the experimentally evaluated laminate properties E^ud = 0.0235 ■ T -f 138.9GPa, 
E,ud = 0.00917 ■ T + 3.93GPa and = 166.17 ■ + 0.206 for T > 400K, 

measured at room temperature after incremental pyrolysis as well as the values for 
GyzUD = 3.36GPa, G^zud = 6.25GPa, Vy:zUD = 0.38, Gmc = 55.7J/m^, evaluated 
at room temperature. 
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Figure 11. a) Segment thickness/width ratio for saturation and b) initiation temperatures 
for segmentation and delamination as a function of the ply thickness 
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The behavior for segmentation initiation and eraek saturation shown in Fig. [TT] was 
observed in the numerie ealeulations. For very thin plies, delamination does not 
oeeur within the temperature range. Transferred to laminate of fabrie this means, 
that the onset of miero delaminations preferentially oeeur at the thieker regions and 
are stopped at the tip of the eross-seetioned fiber bundle. This effeet explains why 
miero delaminations normally do not extend to maero delaminations. 

For a eomparison of experimental with numerie and analytie results, we use a di¬ 
mensionless strain £ = £o/^c = £b, resin = 1.0176 in order to fit the data to the 

temperature eurve. Craeks were eounted only if they span the entire ply and if they 
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Figure 12. Comparison of numerical and analytical approaches with experimental data. 

are elearly visible, using a mieroseope. Therefore, the measured eraek density of 
5.98 1/mm (equivalent to Ngeg ■ t/l = 0.68) has to be eonsidered as a lower bound 
value. Note that the initiation temperatures are only from two experiments and show 
a large seatter. The analytie as well as the numerie approaeh qualitatively eapture 
the fast inerease of the dimensionless segmentation eraek density or segment ratio 
and the saturation in Fig. Reasonable agreement is aehieved between the nu¬ 
merie approaeh and the eraek saturation value of the experiments (ny= 10-20). The 
analytie solution for the saturation density is around 0.45 and far from the experi¬ 
mental value of 0.68. This is mainly due to the model assumption Gmc = Gdc and 
experimental problems of evaluating these values. Nevertheless, the analytie solu¬ 
tion provides a qualitative insight on the role of ply thiekness on initiation values 
(eomp. Fig.fTT]). Note that the segment width in the analytieal model is defined as 
2a — di — c/ 2 . The assumption of nearly eonstant material properties and no ereep is 
believed to be responsible for the differenee in the initiation temperature from the 
experiments. 
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5 Conclusions 


We introduced a disordered spring network model to study the transverse cracking 
of the 90° ply in [0/90] 5 cross-ply laminates, applied for the case of thermal degra¬ 
dation of the matrix material during the pyrolysis process. The main advantage of 
our modeling is that it naturally accounts for the complicated local stress fields 
formed around failed regions in the material, furthermore, it captures the grad¬ 
ual activation of the relevant failure mechanisms and their interactions during the 
fracture process. We have demonstrated that our discrete element model provides 
an insight into the damage process occurring under gradual matrix degradation of 
cross ply laminates. Quantitative results have been obtained on the micro-structure 
of damage and on the development of the porosity for varying ply thickness and 
disorder. 


The results obtained by numerical simulations have been confronted with the ex¬ 
perimental findings and also with an analytic approach. Reasonable agreement was 
found between the numerical, analytical and experimental results. However, the 
numerical simulations proved to be more realistic than the simple analytic ap¬ 
proaches, due to the more detailed description of micro-structure in the discrete 
element model. 


The results of the simulation support material optimization efforts. It was found 
that with decreasing ply thickness, the tendency for the development of delami¬ 
nations is clearly reduced (see Fig. [Tl^). Applied to CMC derived from woven 
fabrics, this could be achieved by reducing the number of filaments per roving, 
the so called K-number. This result gains importance as the development of de¬ 
laminations is the size limiting factor in pyrolytically produced CMC. Thus, larger 
components may be produced by using finer woven laminates. On the other hand, 
this leads to a higher open porosity (see Fig.|^ and in the case of the LSI-process 
to a higher amount of fiber degradation during siliconization. The increase in open 
porosity with decreasing ply thickness can be utilized in producing functionally 
graded CMCs by the LSI-process. A single IK-laminate on top of the standard 
C/C-SiC would give a gradual transition for the SiC-content from the inner parts to 
the totally converted surface of the material. This could give a better adherence of 
coatings like environmental barrier coatings to the substrate. Based on predictions 
for the development of crack morphology, porosity and delamination behavior for 
different system sizes and disorder, promising approaches for new materials have 
been proposed. 
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